Dehydration and ionic conductance quantization in nanopores 
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There has been tremendous experimental progress in the last decade in identifying the structure 
and function of biological pores (ion channels) and fabricating synthetic pores. Despite this progress, 
many questions still remain about the mechanisms and universal features of ionic transport in these 
systems. In this paper, we examine the use of nanopores to probe ion transport and to construct 
functional nanoscale devices. Specifically, we focus on the newly predicted phenomenon of quantized 
ionic conductance in nanopores as a function of the effective pore radius - a prediction that yields a 
particularly transparent way to probe the contribution of dehydration to ionic transport. We study 
the role of ionic species in the formation of hydration layers inside and outside of pores. We find 
that the ion type plays only a minor role in the radial positions of the predicted steps in the ion 
conductance. However, ions with higher valency form stronger hydration shells, and thus, provide 
even more pronounced, and therefore, more easily detected, drops in the ionic current. Measuring 
this phenomenon directly, or from the resulting noise, with synthetic nanopores would provide 
evidence of the deviation from macroscopic (continuum) dielectric behavior due to microscopic 
features at the nanoscale and may shed light on the behavior of ions in more complex biological 
channels. 
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I. INTRODUCTION 

The behavior of water and ions confined in nanoscale 
geometries is of tremendous scientific interest. On the 
one hand, biological ion channels, which form from mem- 
brane proteins, perform crucial functions in the cell 
[U H]. On the other hand, there have been recent ad- 
vances in aqueous nanotechnology such as nanopores 
and nanochannels, which hold great promise as the 
basic building blocks of molecular sensors, ultra- fast 
DNA sequencers, and probes of physical processes at 
the nanoscale [3]. Indeed, nanopore-based proposals for 
DNA sequencing range from measuring transverse elec- 
tronic currents driven across DNA [4-_8^ to voltage fluc- 
tuations of a capacitor [9HTT] to ionic currents [I2HI9] . 

Recent experiments show that we are tantalizingly 
close to realizing a device capable of ultra- fast, single- 
molecule DNA sequencing with nanopores: identifica- 
tion of individual nucleotides using transverse electronic 
transport [20, 2Ij has been demonstrated. Discrimination 
of nucleotides using their ionic blockade current when 
driving them individually though a modified biological 
pore has also been demonstrated [22l [23] . In these sys- 
tems, the presence of water and ions will affect the signals 
and noise measured and thus understanding their behav- 
ior is an important issue in both science and technology. 

Many computational studies have been dedicated to 
relating the three-dimensional structure [I] |211 [55] of bi- 
ological ion channels to their physiological function, e.g., 
ion selectivity. For instance, recent studies have exam- 
ined the role of ligand coordination in potassium selec- 
tive ion channels [25^29] . Biological channels, however, 
are complex pores with many potential factors contribut- 
ing to their operation. Thus, only in a limited number 
of cases have universal mechanisms of ion transport been 
investigated, such as the recent work on the role of "topo- 



logical constraints" in ligand coordination 




Figure 1: Schematic of ion transport in the presence of hydra- 
tion layers. Only the first hydration layer is shown for sim- 
plicity. Ions in bulk water form hydration layers that make 
the ion behave as a "quasi-particle" that includes the ion and 
tightly bound water molecules. This quasi-particle is then 
solvated in the high-dielectric water. As the ion goes from 
the bulk solution to the pore it has to partially shed its hy- 
dration layers, i.e., the quasi-particle has to break apart. This 
gives a nonlinearity in the energetic barrier to transport. A 
continuum picture neglects these features and considers only 
the dielectric barrier that ions have to overcome by moving 
from bulk water with dielectric constant e^ into the inho- 
mogeneous, low-dielectric pore environment with tp ^ e™. 
Specialized proteins facilitate this process in biological ion 
channels via the presence of localized surface charges/dipoles 
and other mechanisms. 

Fundamental developments in the fabrication of syn- 
thetic nanopores j53 40 , however, open new venues for 
investigating the behavior of ion channels and dynami- 
cal phenomena of ions, (bio-)molecules, and water at the 



nanoscale. For instance, what are the dominant mecha- 
nisms determining ionic currents and selectivity? What 
role do binding sites play versus hydration in constrained 
geometries? How accurate are "equilibrium" and/or con- 
tinuum theories of ion transport? Well-controlled syn- 
thetic pores can be used in this context to examine how 
ion transport is affected, for instance, by changing only 
the pore radius, in the absence of binding sites and sig- 
nificant surface charge within the pore. 

In this paper we examine the role of dehydration in 
ionic transport through nanopores. In particular, we in- 
vestigate the recent prediction of quantized ionic con- 
ductance by two of the present authors (MZ and MD) 
pi] . namely that drops in the conductance, as a func- 
tion of the effective pore radius, should occur when suc- 
cessive hydration layers are prevented from entering the 
pore. This effect is a classical counterpart of the elec- 
tronic quantized conductance one observes in quantum 
point contacts as a function of their cross section (see, 
e.g., Ref. [42j). We examine different ions, both positive 
and negative, and of different valency (namely, Cl~ , Na"*" , 
K+, Ca^+, and Mg^+). We find that the ion type plays 
only a minor role in the radii of the hydration layers, 
and thus does not affect much the pore radii at which a 
sudden drop in the current is expected. Divalent ions, 
however, are the most ideal experimental candidates for 
observing quantized ionic conductance because of their 
more strongly bound hydration layers. Further, the fluc- 
tuating hydration layer structure and changing contents 
of the pore should give a peak (versus the effective pore 
radius) in the relative current noise - giving an addi- 
tional method to observe the effect of the hydration lay- 
ers. Thus, we elucidate how quantized ionic conductance 
provides a novel tool to deconstruct the energetic contri- 
butions to ion transport. 

The paper is organized as follows: In Sec. |llj we give 
a macroscopic (i.e., a continuum electrostatic) viewpoint 
on the energetics of ion transport. In Sec. |III[ we exam- 
ine how ions induce local structures in the surrounding 
water known as hydration layers - an effect that is not 
taken into account when using continuum electrostatics 
to estimate energetic barriers to transport. Further, we 
calculate the energies stored in these layers and develop 
a model for the energetic barrier for ions entering a pore. 
In Sec. |IV[ we use a Nernst-Planck approach to relate 
this barrier to the ionic current. In Sec. |V] we discuss 
how the presence of the hydration layers gives rise to a 
peak in the relative noise in the ionic current at values of 
the effective pore radius congruent with a layer radius. 
In Sec. |VI[ we then present our conclusions. 



II. IONIC TRANSPORT 

The experimental set-up we are interested in is that 
of ions driven through a pore/channel of nanoscale di- 
mensions under the action of a static electric field [71] . 
Such a situation is depicted in Fig. [l] A simple ap- 



proach to ionic transport is to envision the ions moving 
through an energetic barrier due to going from the high- 
dielectric aqueous environment into the inhomogeneous, 
low-dielectric environment of the pore, treating the sur- 
roundings as continuum media. The resulting approach 
is inherently static: by analyzing the energetic barrier 
to (near-equilibrium) transport one obtains information 
about how different factors - the pore material (through 
its dielectric constant), the pore dimensions, the presence 
of surface charges, and the presence of the high-dielectric 
water along the pore axis - would affect transport. 

Indeed, one of the first calculations of the dielectric 
barrier (using a "Born solvation" model) was done by 
considering the ion solvated in water and moved into a 
low-dielectric, pore-less membrane |33J|33]. This provides 
an estimate of the energies involved by calculating the 
energy change of solvating the ion in continuum water, 
with dielectric constant e^ ~ 80, to "solvating" it in a 
continuum material with e^ ~ 2 (representative of lipid 
membranes [75]). For instance, the energy change of a 
Cl~ ion, with effective radius R^2A [73] , moved from 
continuum water to the continuum material is 
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This is quite a substantial energy change - about half the 
solvation free energy of Cl~ [H US] . The finite thickness 
of the membrane does not change this value significantly. 
For thick membranes, it is lowered by [JS] |31] 
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for ew 3> Ep, where I is the membrane thickness (and pore 
length). For e^ ~ 2 and e„ « 80, this gives ~ 5/1 eV A « 
0.1 eV for a membrane of thickness I — 50 A. That is, 
the Born estimate in Eq. (fTl) is lowered to ^ 1.7eV. 
However, the membrane width [46^ and composition can 
play a significant role in this estimate. For the common 
synthetic pores made of silicon dioxide (cp w 4) or silicon 
nitride {cp ~ 7.5) the estimate in Eq. M is reduced from 
~ 1.8 eV to ^ 0.9 eV and ~ 0.4 eV, respectively. These 
barriers are more than an order of magnitude larger than 
k^T at room temperature, where kg is the Boltzmann 
constant. 

Due to this magnitude, it is clear that the energy scale 
of solvation is one of the controlling factors in ion trans- 
port. However, in addition to the above there is water 
present in the pore. One expects, therefore, that the 
energy of solvation would be decreased from simple es- 
timates like that of Eq. Q. Several groups have cal- 
culated this contribution pBl H^ HBVH5] . For instance, 
Ref. [in] shows that the energy barrier of bringing an 
ion from continuum water into a low-dielectric, contin- 
uum membrane is reduced from ~ 40kcal/mol ~ 1.7 eV 
to ~ 20kcal/mol w 0.9 eV by the presence of water in 
the pore. This demonstrates that a pore filled with a 



high dielectric medium (e.g., continuum water) can sig- 
nificantly lower the barrier to transport. Even still, the 
barrier remains substantial. 

In biological systems, however, the pores provide a 
channel with a much lower barrier as indicated by the 
conductance of many biological ion channels. These 
pores are formed from specialized proteins whose role is 
precisely to facilitate passage of ions (and further to selec- 
tively allow passage of certain ions). Clearly, pores with 
internal charges and/or dipoles can significantly reduce 
the energetic barriers for transport. Indeed, the effect of 
surface charges has been calculated in clean pores [50l - 
152] and when present in sufficient amounts would negate 
the effect we predict as the reduction of the energetic 
barrier would be comparable to, or larger than, the hy- 
dration layer energies. Therefore, our interest is in clean 
pores with little to no surface charge where clear-cut ex- 
periments can be performed to understand the effect of 
hydration on transport. This rules out the direct use of 
some biological ion channels, particularly those with very 
small pores where single-file transport occurs [211 [531 IM] , 
because of the presence of localized charges and dipoles. 

To conclude this section, we note that the contin- 
uum description suffers from a number of issues at the 
nanoscalc: it is only valid beyond the correlation length 
of the material |55| , which for the strong fields around an 
ion is ~ 8 A for water (see below), similar to the ~ 5 — 8 A 
in water only |56j ; linear continuum electrostatics is only 
valid when the polarization field is co-linear with the elec- 
tric field (not the case in the hydration layers we dis- 
cuss below); in a related issue, it is only valid for weak 
fields (in the context of ion channels, see, for example. 
Sec. 3.4 in Rcf. [49J); there is also an issue of where the 
"surface" separating the charge and the dielectric mem- 
brane/continuum water is located, especially for fiuctuat- 
ing atomic ensembles as is the case for protein pores and 
molecular (rather than continuum) water. Thus, while 
a continuum picture can highlight some general features 
of the energetic barrier to ion transport - in some cases 
giving compact analytical expressions - it breaks down 
when trying to understand the effect of structure at the 
nanoscalc. In fact, macroscopic, continuum electrostatics 
is not designed to study specific features or short-range 
interactions at these length scales. This is precisely what 
we seek to address in the following sections. 



III. HYDRATION OF IONS 

We begin our study of quantized conductance by first 
illustrating how ions are hydrated in solution and then 
discuss the energies involved in this process. The forma- 
tion of hydration layers around ions has been known for 
some time (see, e.g., Ref. [I]), and is due to the strong 
local electric field around the ion and to repulsive short- 
range interactions among molecular/atomic species. We 
use molecular dynamics (NAMD2 [57j) simulations to un- 
derstand the structure of hydration layers when different 



ions are inside and outside of nanopores [71] . 

Figures [2] and [3| show the water density oscillations for 
several common ionic species |75) . There is a strong peak 
in water density about 3 A away from the ions, with two 
further oscillations after that spaced about 2 A apart. 
These oscillations signify that there are strongly bound 
water molecules forming around the ions. Table [ij lists 
the hydration layer radii from both this study and ex- 
periment. We find very good agreement with the experi- 
mental data for all cases. The water density approaches 
the bulk value (- 0.033 Molecules/A^) at about 10 A, 
which is also consistent with the experimental value. 

The oscillations in water density also give rise to oscil- 
lations in the local electric field. Figures[2|and[3|show this 
for monovalent and divalent ions where the time-averaged 
electric field was calculated from the bare ion value plus 
a sum over all partial charges given by the hydrogen and 
oxygen atoms of water ^6]. In the figures, the first hy- 
dration layer gives pronounced field oscillations for all 
species examined. The other oscillations in the field are 
more well-defined for K+, Ca^^, Mg^^, and to some ex- 
tent Na^, compared to Cl~. Anions, such as Cl~, have 
a different structure of the water around them compared 
to cations: in the first layer, they pull one of the hydro- 
gen atoms of each of the water molecules closer while the 
other interferes with the formation of the second layer, 
possibly hindering the ability of the second layer to form 
a "perfect" screening surface. The fact that the electric 
field is not simply suppressed by l/e^u shows the diffi- 
culty of a macroscopic (continuum) dielectric picture to 
predict behavior at the nanoscale (similar to well-known 
features in other systems such as Friedel oscillations and 
apparent from the derivation of continuum electrostat- 
ics, where averaging is required over length scales much 
larger than the correlation length of the material [55]). 

We now estimate the energies contained in these layers, 
which we list in Table |l[ The electric fields seen in Figs. [2] 
and [3] show an oscillating behavior that is reminiscent of 
a set of Gauss surfaces, i.e., layers of alternating charge 
that screen the field of the ion. Thus, in order to esti- 
mate the energies contained in the layers, we replace the 
microscopic structure giving rise to the complex field by 
a set of surfaces as shown in Fig. [4] that perfectly screen 
(with dielectric constant e^,), rather than over-screen, the 
ion charge. 

Within this picture, the energy of the i*'' hydration 
layer of ionic species ly is [H] 
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where qi, is the ionic charge and R^^ are the inner 
(outer) radii demarcating the hydration layer as obtained 
from the water density oscillations. In order to obtain the 
innermost radius we compute the total solvation energy, 
Ut — —qt [f-w — Ep) / (STreoepCiu-Ri^), and compare with 
the experimental free energies [45] . which are dominated 
by the electrostatic energy. These free energies, together 
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Figure 2; Top panels: Water density oscillations versus distance for Cl~, Na"*", and K"*" in bulk water. Black, solid lines indicate 
the density calculated from the oxygen atom positions for Cl~ and hydrogen atom positions for the cations. The arrows 
indicate the minimum in the density oscillations. The blue, dashed lines indicate the density calculated from the hydrogen 
atom positions for Cl~ and oxygen atom positions for the cations. Bottom panels: The electric field due to both the bare ion 
(red, dashed line) and due to the ion plus partial charges on the water molecules (black, solid). The arrows again indicate the 
minimum in the density oscillations. 
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Table I: Table of physical quantities from simulation and experiment. The theoretical hydration layer radii, Ri, for all three 
layers are defined using the ion-oxygen distance for Cl~ and ion-hydrogen distance for the cations. The first oxygen density 
maximum, Rxo, is for all ions X using the present theory and experiment (average values taken from Ref. |3S]). The second 
maximum is also shown for the divalent ions from both theory and experiment. The inner/outer radii that enter Eq. Q are 
shown, the first of which is calculated such that Eq. (J4} equals AG (exp) when _R^ — >■ oo (see also text). The next three 
inner/outer radii are taken from the minima of the oxygen density for Cl~ and the minima of the hydrogen density for the 
cations. Further, we report the layer energies Ui (using ep = 1), and the Gibb's free energy from experiment [45], and the 
experimental mobilities ITI used in this work. 



with the layer energies (for Cp — 1), are tabulated in Ta- 
ble |Tj Except for the third hydration layer for monovalent 
ions, the layer energies are greater than other free energy 
contributions such as the entropy change due to the water 
structure or van der Waals interactions [SHI [50] . In Eq. 
^ we have also added a possible screening contribution, 
Cp, from the pore material and/or charges on the surface 
of the pore. In Ref. ^T] this was assumed to be one: 
the low-dielectric pores reduce the energy barrier only 



by a small amount and in a different functional form. In 
Sec. |V]we will discuss the effect of this screening on the 
detection of quantization steps. 

Previously, we proposed a model for how the energy 
is depleted in a hydration layer as the effective radius of 
the pore, Rp, is reduced [41]. In this model, the energy 
change is proportional to the remaining surface area of 
a hydration layer within a pore. It takes into account 
both that the water-ion interaction energy of small wa- 
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Figure 3: Top panels: Water density oscillations versus distance for Ca^""" and Mg^**" in bulk water. Black, solid lines indicate 
the density calculated from the hydrogen atom positions. The arrows indicate the minimum in the density oscillations. The 
blue, dashed lines indicate the density calculated from the oxygen atom positions. Bottom panels: The electric field due to 
both the bare ion (red, dashed line) and due to the ion plus partial charges on the water molecules (black, solid). The arrows 
again indicate the minimum in the density oscillations. 
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ical shells when an ion is inside a pore of small enough 
radius (see Fig. 1 in Ref. [H]). Contributions from, e.g., 
van der Waals interactions with the pore and changes in 
the water- water interaction, are small [SHIBD]. Thus, the 
energy of the remaining fraction /^i, of the i*'' layer in 
the pore is taken as JJw = fwU°^. The fraction of the 
layer intact is fi^, = Sii^/AirRf^ where Si^ is the surface 
area (of the spherical layer) remaining where the water 
dipoles can fluctuate. The latter is given by 
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dORi 
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where O (x) 
sin~^(i?p/i?i^). 
face left is 



is the step 
When Rp < R., 



function and 9ci, = 
I, , the fraction of the sur- 



Figure 4: The magnitude of the approximate electric field 
(given by a set of Gauss surfaces, see inset) used to represent 
the fields in Figs. [2] and [S] and also used to calculate the 
energy contained in each layer (Eq. Q). 



ter clusters is approximately linear in the number of wa- 
ters [601 I61j and that molecular dynamics simulations 
show a time-averaged water density with partial spher- 
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The total internal energy change will then result from 
summing this fractional contribution over the layers to 
get 
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Sec. IV below), we calculate the free energy [7S] change 



for species v as 



AF, = AC/, - TA5,, 



(8) 



which includes an entropic contribution from removing a 
single ion from solution and localizing it in the pore re- 
gion. This entropic contribution is A5,y = fcBln(P^rio), 
where we have assumed an ideal ionic solution and Vp is 
the volume of the pore and no is the bulk salt concentra- 
tion for all species v. The free energy change is plotted 
in Figs. [6] and [7] versus the effective pore radius and it 
is substantial when the latter becomes smaller than the 
outer hydration layer. 



Figure 5: Distributions of the dipole orientations of water 
molecules within a cylindrical annulus 5 A away from, i.e., 
in the second hydration layer of, a Ca^"'' ion in bulk (black 
line) and inside pores of radius 8 A (red dashed line), 12 A 
(green dash-dotted line) and 15 A (blue dotted line). The 
mean value is around 0.38 (corresponding to the water dipole 
pointing 68 degrees away from the ion- water vector), except 
for the 8 A pore, which increases to 0.42 (corresponding to 
the water dipole pointing 65 degrees away from the ion- water 
vector). This signifies a moderate tightening of the water 
dipole around the field of the ion as the pore size is reduced. 



We stress first that the effective radius R^^ is not nec- 
essarily the nominal radius defined by the pore atoms. 
Rather, it is the one that forces the hydration layer to be 
partially broken because it can not fit within the pore, 
and it could be smaller than the nominal pore radius 
by the presence of, e.g., a layer of tightly bound water 
molecules on the interior surface of the pore. Second, 
our model misses internal features of the hydration lay- 
ers themselves. For instance, Ref. i62j examines the first 
hydration layer structure in carbon nanotubes of different 
radii. These authors find a large increase in the energy 
barrier when the pore radius nears the inner hydration 
layer. They also seem to observe sub-steps in the wa- 
ter coordination number within the inner shell as the 
pore radius is reduced. Thus, although our model con- 
tains only a single "smoother" step, experiments could 
very well observe these internal sub-steps corresponding 
to the sudden loss of a single or few water molecules out 
of a given hydration layer. 

Another basic assumption in our model is that the in- 
teraction energy of the water molecules in each layer is 
the same regardless of whether the ion is inside or outside 
of the pore. Figure [5] shows the distribution of the dipole 
orientation of water molecules both in bulk and inside 
pores of different radius [77] . The average dipole orienta- 
tion of the waters changes very little inside the pore, as 
do their fluctuations, thus supporting this assumption. 
In addition, the structure of the first two hydration lay- 
ers (not shown) is essentially the same in and out of the 
smallest (8 A) pore. 

In order to make a connection with the ionic current (in 



IV. IONIC CURRENTS 

We now want to relate these energy barriers to the 
ionic current through the pore [75]. We do this by solv- 
ing the Nernst-Planck equation in one dimension. Since 
this model consistently solves for both drift and diffu- 
sion contributions to ionic transport, and yields a com- 
pact analytical expression, we use it below with the en- 
ergetic barriers found from the above model of dehydra- 
tion. Even though this analytical model does not include 
some effects such as ion-ion interaction, we expect that 
it is qualitatively accurate as discussed along with its 
derivation. 

The steady-state Nernst-Planck equation (see, e.g., 
[H 1631164]) for species u in one dimension (assuming vari- 
ability on the pore cross-section is not important) is 



Jiy = -qvD^ 



dn^ (z) 
dz 
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dz 
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where Ji, is the charge flux for species v, z £ (0, 1) 
is the axial coordinate along the pore axis of length 
I, n^ [z) is the ion density, Di, is the diffusion coefh- 
cient (assumed to be position independent), and $, (z) 
is the position-dependent potential (including both elec- 
trostatic and other interactions that change the energy 
within the pore). A full solution would require solving 
the density and potential within the reservoirs and pore 
simultaneously (see, e.g., Ref. [65]). However, we deal 
with high-resistance pores. Thus, we approximate the 
left [L) and right [R) reservoirs with constant concen- 
trations riL and n/j, and the boundary conditions at the 
edge of the pore are rii, (0) = n^ and rii, (Z) = nji. This 
is equivalent to assuming that as soon as an ion leaves or 
enters the pore, the ions in the immediate surroundings 
of the pore equilibrate rapidly to their prior distributions. 
Thus, multiplying by exp [q^^u {z) /kBT) to get 

dz L J 

and integrating yields the flux for species i^ as 



Ji> = -qwDi 



njie 



qu<i>u{l)/kBT 



riLe 



g,.*.(0)/fci5T 



/^ dze«^*^ 



{z)/kBT 



(11) 
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Figure 6: Free energy changes, AF^, and currents versus the effective pore radius for the monovalent ions and a field of 
i mV/A. The black, dashed line indicates the free energy change. The remaining lines indicate the current with different 
standard deviations of the noise (see text for details). The currents are for a = P^F, with p = 0.3, 0.2, 0.1, from top to 
bottom. 
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The black, dashed line indicates the free energy change. The remaining lines indicate the current with different standard 
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We make the further simplifying assumption that the 
electrostatic potential drops linearly over the pore - rec- 
ognizing that in the presence of a significant potential 
barrier, e.g., due to the stripping of the water molecules 
from the hydration layers and in the absence of sur- 
face/fixed charges in the pore, the ionic density in the 
pore is small and thus the field is due to ionic charge lay- 
ers on both sides of the pore. Results from many works 
that include ion-ion interactions indeed find a linear drop 
of the potential across the pore (see, e.g., Ref. [66]). In 
this case, ions form a capacitor across the pore and ev- 
ery so often one ion translocates through the pore. The 
"heaUng" time for the loss of this ion is very short [66] 
and, thus, the field (potential drop) is not strongly af- 
fected [HO]. Also, we assume that the potential barrier 
due to changes in these other interactions is constant over 
the pore - this ignores a region near the pore entrance, 
but will not qualitatively change the solution. Therefore, 
the potential for species v can be written as 



<J'.W 



V . AF, 



(12) 



when z € (0, 1). The boundaries are given by $2 (0) — 
and $2 (l) = V . Performing the remaining integral and 
for equal reservoir densities (our case), nj^ — nji — uq, 
we get 



J.= 



IksT 



(13) 



Relating the diffusion coefficient to the mobility via the 
Einstein relation, /i^ = q^D^/kBT , and putting in the 
constant electric field E — V/l, one obtains 



Ju - -q.nofi.Ee'^^"/''^^. 



(14) 



That is, the flux of an ionic species is proportional to the 
electric field and density, where the latter is suppressed 
by a Boltzmann factor [H] . 

Now that we have an expression relating the energy 
barrier to the transport properties, we can calculate the 
current as a function of effective pore radius by multiply- 
ing Eq. ( 14 ) by the area of the pore to get 



2Tr R;J^ 



'i/oe 



^AF^/ksT 



(15) 



where we have defined a standard current I^q — 
—q,yno2TrRpiJ,yE that would flow in the absence of an en- 
ergy barrier. The current (151, with the mobihties and 
energies in Table III along with Eqs. Q-lls]), is plotted 
in Figs. |6] and [7] as a function of effective pore radius 
and for a field of ImV/A [81 . The energetic barriers 
create sudden drops when the pore radii are congruent 
with a hydration layer radius. These correspond to the 
quantized steps in the conductance. 



V. EFFECT OF NOISE 

In a real experiment, there will also be fluctuations in 
the energetic barrier due to the fact that the hydration 
layers are not defined by their time-averaged value (i.e., 
they are not perfect spherical shells) and also due to fluc- 
tuations of the water structure and contents of the pore 
(both within a single experiment and also structural vari- 
ations between experiments). Thus, we also examine the 
effect of these fluctuations and the current noise they in- 
duce. Thus, we calculate an averaged current for species 
V as 



where the factor A is 



(I. 



{hoe 



-AF^/kBT\ 



(16) 



We consider two speciflc models: Gaussian fluctuations 
of the free energy with a standard deviation proportional 
to the free energy barrier at a flxed pore radius and 
Gaussian fluctuations in the effective pore radius. The 
latter was also considered previously [H] where it was 
found that this type of noise smooths out the visibility of 
the drops in conductance (i.e., the peaks in the deriva- 
tive d{Ii,)/dRp become smoother with increasing noise). 
However, it was also shown that this fluctuation induces 
a peak in the relative current noise that is much less sen- 
sitive to the strength of the fluctuations - thus giving an 
alternative method to detect the effect of the hydration 
layers. We develop a model for this relative noise here 
but do not perform the calculation of Eq. ( 16 ) for all the 
different species. 

Fluctuating energy harrier - The first model we con- 
sider is an energy barrier that fluctuates according to a 
Gaussian distribution. We neglect fluctuations that make 
the barrier negative, so that the average current is 



(/,) - ^ / d{AF) e-AFA-BTg-(Af-AF.)V2a^ 
■^o- Jo 



(17) 
where a is the standard deviation of the fluctuations and 

/•OO 

Jo 
is the normalization. The average current is thus 



(/,)=/.oAe-(^^--'^'/2'=«^)/'=-^, 



(19) 



erfc 



A^ 



(-AF, + aVfcsT) /\/2^ 



erfc 



(20) 



The value of A for small a is very close to 1. Thus, the 
effect of a fluctuating energy barrier with small fluctua- 
tions is simply to lower the energy barrier by an amount 
a^ /2kBT. For stronger fluctuations, the factor contain- 
ing the complementary error function, erfc, gives different 
limiting dependencies of the average current as the fluc- 
tuation strength a is increased. However, large fluctua- 
tions are well outside the realm of validity of the present 
model. 

The relative noise in the current provides even more 
information. The relative noise is 



A/,, 



{I) 



(21) 



The expectation value of the square of the current is given 
by 



T fOO 

(/2) ^ ^ / rf(AF) e-^^/fc-^e-(^^-^^'^)'/^^' 

•^(T Jo 



= BI^ g-(2AF„-2<TVfcBT)/fcBT^ 



(22) 



Where the normalization is as before and the factor B is 
given by 



erfc 



B 



(-AF,, + 2aykBT) /V2a^ 



erfc ~AF^/V2^ 



(23) 



Thus, the relative current noise induced by an energy 
barrier with fluctuations is 



A/,e, = A/e-V('c..T)^^_l. 



(24) 



For small fluctuations, A and B depend very weakly on 
a and are both very close to 1, giving a relative current 
noise 



Alrel ~ a/ksT. 



(25) 



As expected, the relative noise increases with the 
strength of the fluctuations. For fluctuations propor- 
tional to the energy barrier, as shown in the Figs. IHland 
[Tj the fluctuations give rise to a monotonic increase in the 
relative noise. Overall, the effect of fluctuations in the en- 
ergy barrier is to decrease the effective energy barrier and 
increase the current. This reduces the magnitude of the 
drops in the conductance but does not destroy their vis- 
ibility. This would therefore help in observing quantized 
ionic conductance. It is worth noting, however, that this 
type of noise makes the step of the third hydration layer 
the most pronounced. This seems an unlikely situation 



in actual experiments and other types of noise need to 
be considered. 

Fluctuating effective pore radius - In addition to the 
above noise, one expects that there would be fluctua- 
tions in the radii of the hydration layer/nanopore system. 
Previously, we demonstrated that this type of noise can 
smear the effect of the steps in the current [JT] . As was 
seen, however, this noise also gives a peak in the rela- 
tive noise in the current that is much less sensitive to the 
fluctuations than the average current. Here we develop 
a model of this behavior by calculating the relative noise 
assuming fluctuations across a single, perfect step in the 
free energy (see the inset of Fig. [s]). 

The average current due to species v when averaged 
over fluctuations in the effective pore radius is 



The relative noise is 



1 



dRI^oiR)e 



-AF(fl)/fcBTg-(fl,-flp)V2?" 



(26) 

where ^ is the standard deviation of the radial fluctua- 
tions, A/j is the normalization, and the explicit R depen- 
dence has been included in both the barrier AF and the 
prefactor I^q. The dominant factor is the exponential 
of the free energy barrier and the quadratic dependence 
of /yo on R can be ignored. For small fluctuations, the 
lower limit of the integral can be extended to — cx) and 
Af^ -^ \J2'kC to give 



luojRp) 



dRe 



~AF{R)/kBT^-{R-R„y/2i^ 



(27) 

Previously, we performed the averaging according to Eq. 
(26) [41^, but here we instead use Eq. (271 with the 



approximate energy barrier AF (R) = AF^Q {Rh — R) 
of a single hydration layer of radius Rt and take I^q to 
be the current in the absence of the barrier. The average 
current then becomes 



(/.) = I 



uO 



^-AF^/kBT 



il-C) + C 



where 



C 



Irfc f^^L^ 
2 V V2^ 



Similarly, for the square of the current one finds 



(28) 



(29) 



(30) 



Although (/j^) and (/^) are dependent on the strength of 
the fluctuations, ^, the relative current noise has a uni- 
versal behavior in the parameter R = [Rh — Rp) /\/2^. 
That is, all features in the relative noise would be present 
regardless of the strength of the noise. However, the peak 
in the noise (see below) shifts to smaller values of Rp as 
the noise strength is increased, which is qualitatively in 
agreement with the full averaging (Eq. (26)) performed 
in Ref. m]. 



^J {l - c-^F'^/'^btY c {1 - C) 

For large or small Rp, the relative noise goes to zero, 
which can be seen from the properties of erfc that make 
C — >■ 1 and C — >■ for large and small Rp, respectively. 
In between these limits, there would be nonzero relative 
noise, therefore indicating that the relative noise would 
have a maximum. The peak in the relative noise occurs 
for Rp < Rh- For a large energy barrier AFh, this peak 
occurs when C is small. Thus, we can approximate the 
relative noise as 



A/rcl 



C 



-AFh/ksT _|_ (J- 



This gives a peak in the noise when C : 
a value 



(32) 



.-AF./ksT ^i^^ 



A/,^1 



1 



,AFh/2kBT 



(33) 



The peak is exponentially large in the energy barrier. 
However, the model with the electrostatic energy given 
by Eq. ([7| does not have an ideal step in the free energy 
(see, e.g., Figs.l6]andl7]). From previous work "41], we can 
identify the peaks, R*, and use AFh ~ AF^, (i?*) . This is 
done in Fig. [sjfor Cl~. The model agrees quantitatively 
with the full averaging performed in Ref. [HI. The only 
feature missing is the additional background noise away 
from the step due to the non-uniform energy barrier on 
both sides of the step. 

Thus, from this "two-channel" model of noise we have 
found two generic features: (i) a peak develops in the 
relative current noise that is exponentially high with the 
hydration energy barrier, and (ii) it is present regard- 
less of the noise strength, although its location moves to 
smaller values of the pore radius with increasing noise 
(likewise, the peak becomes wider). These features are 
in agreement with what is found from performing the full 
averaging from Eq. ( 26 ) using the surface area model of 
the energy barrier [41j . In the full model the fluctua- 
tions will eventually smooth out the peak in the relative 
current noise. The latter, however, is still much less sen- 
sitive than the average current drops, making the peak 
in the relative current noise versus Rp a robust indicator 
of dehydration. 

Barrier reduction - In addition to the above fluctua- 
tions, there are factors that reduce the energetic barrier, 
such as the presence of some surface charge and/or dielec- 
tric screening in the pore. In Eq. (H) we included a dielec- 
tric constant ep to represent a reduction in the hydration 
layer energy barrier from these sources. We expect, how- 
ever, that the introduction of this constant overestimates 
the barrier reduction. It amounts to replacing the water 
molecules screening the ion with a material of lower di- 
electric constant but in the exact geometry of the water 
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factors under control greatly affects the transport prop- 
erties of the ions [i67j. Therefore, pores made of, for in- 
stance, semiconducting nanotubes may be ideal. Indeed, 
pores made of these materials have been recently demon- 
strated [6S]. However, rough surfaces that are present 
in pores made of, e.g., silicon nitride, should still allow 
for quantized conductance to be observed, so long as the 
variation of the effective radius of the pore is not too 
strong. The noise in the effective radius of the pore was 
investigated previously in Ref . [41] , where we found that 
only beyond variation in the radius of 0.2 — 0.3 A will 
the effect be washed out. However, even beyond this 
variation magnitude, the relative current noise signifies 
the presence of steps in the energy barrier, thus giving a 
more robust indicator of the hydration layers' effect on 
transport. 



Figure 8: Relative current noise induced by structural fluc- 
tuations in the efi'ective pore radius and/or hydration layer 
radius. The inset shows the approximate change in free en- 
ergy A_F)i as a function of effective pore radius in proximity 
to a hydration layer structure. The noise in the pore ra- 
dius induces fluctuations between the high and low energy 
states. Here the third hydration layer radius of Cl~ is taken, 
i?h — 7.1 A. See the text for details on A_Fii. The fluctuation 
strength from right to left is ^ = 0.05, 0.15, 0.25, 0.35, 0.45 A. 



molecules. This is very unlikely since the pore screening 
comes from the fixed surface of the pore and thus in a 
different functional form. Nevertheless, it is instructive 
to see how the drops in the current are reduced by this 
effective lowering of the energy barrier. Figures [9] and [lO] 
show the energy barrier and current for several values of 
this effective dielectric constant. We find that even for 
fairly large Cp (~ 7) , the barriers are large enough to give 
a noticeable drop in the current. 

Bulk concentration - We also mention the effect of 
changing the concentration of ions in bulk. We have as- 
sumed that the hydration layers are well formed away 
from the pore. Large ionic concentration in bulk, how- 
ever, would affect the formation of the hydration layers. 
For a completely disassociated 1:1 salt, the ion- ion dis- 
tance goes as ^ S^AJUq A where the bulk concentration, 
no, is given in mols/L. Thus, the inter-ion distance is 
~ 9.4 A for a 1 M solution, which is almost large enough 
to house both the first and second hydration layers. How- 
ever, concentrations lower than 1 M are preferable. 

Some remarks - We have discussed many of the factors 
that will affect the detection of quantized ionic conduc- 
tance. The most ideal experiment would be one with 
pores of well-controlled diameter and with smooth sur- 
faces. Likewise, a small (or no) amount of surface charge 
and a low dielectric constant of the pore will make the 
effect more pronounced (and the ability to gate a pore, 
e.g., made of a nanotube, would help even more in under- 
standing the energetics of transport). Not having these 



VI. CONCLUSIONS 



Ionic transport in nanopores is a fascinating subject 
with a long history and impact in many areas of science 
and technology. Recent work on developing aqueous- 
based nanotechnology and understanding biological ion 
channels requires a firm understanding of how water 
and ions behave in confined geometries and under non- 
equilibrium conditions due to applied fields. For exam- 
ple, the quest for ultra- fast, single-molecule DNA se- 
quencing has yielded a number of proposals based on 
nanopores 0. Among them, transverse electronic trans- 
port m [S] (whose theoretical basis includes the investiga- 
tion of atomistic fluctuations [IHZj and electronic noise in 
liquid environments j8j) and ionic blockade p^HT9| have 
yielded promising recent experiments (Refs. HOI HI] and 
[351153], respectively). In all these cases, both water and 
ions are present and will have a signiflcant impact on the 
signals and noise observed. 

In this work, we have analyzed in detail the recent 
prediction of quantized ionic conductance |41| and ex- 
amined how different aspects of the ion-nanopore system 
influence the detection of this phenomenon. Namely, we 
have shown that the ion type affects very little the radii 
at which the conduction should drop. High valency ions, 
however, should give even more pronounced drops in the 
current and thus may help in detecting this effect. Fur- 
ther, the presence of the hydration layers gives a peak 
in the relative noise at pore radii congruent with a layer 
radius. This relative noise is much less sensitive to fluctu- 
ations than the average current, and provides a promising 
approach to detect the effect of hydration. 

Overall, quantized ionic conductance yields experimen- 
tal predictions that will shed light on the contribution of 
dehydration to ion transport and we hope this work will 
motivate experiments in this direction. 
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Figure 9: Free energy changes, AF, and currents versus the eflective pore radius for the monovalent ions and several values of 
Ep (tp = 7, 5, 3, 1 from bottom to top) and for a field of 1 mV/A. The dashed lines indicate the free energy change. The solid 
lines indicate the current (see text for details). The currents are for ep = 7, 5, 3, 1, from top to bottom. 
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Figure 10: Free energy changes, AF , and currents versus the effective pore radius for the divalent ions and several values of 
Ep (tp = 7, 5, 3, 1 from bottom to top) and for a field of 1 mV/A.The dashed lines indicate the free energy change. The solid 
lines indicate the current (see text for details). The currents are for ep = 7, 5, 3, 1, from top to bottom. 
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Although similar conclusions should apply in other sce- 
narios such as the generation of a concentration gradient. 
In pores - especially biological pores - the membrane di- 
electric constant, em, and the dielectric constant of the 
pore material, tp, can be different. 

The effective radius can be estimated from, e.g., molec- 
ular dynamics simulations that give a surface where the 
screening charge due to the hydrogen or oxygen atoms 
of water fluctuates. For instance. Figs. [2] and |3] show this 
surface (see also Refs. [691170] ). 

For an ion in bulk water, we simulated a hexagonal box 
of 150 A height and 43 A radius with periodic boundary 
conditions in all directions. We then fixed an ion in the 
center of the box and counterion(s) near the edge of the 
box, far away from the ion of interest. For an ion in a 
pore, a cylindrical pore of radius R was cut into a hexag- 
onal silicon nitride film 97 A thick and of 29A radius. 
This was accomplished by removing all silicon and nitro- 
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gen atoms within a distance R from the z-axis. An ion 
was fixed in the middle of the pore and counterion(s) were 
fixed outside of the pore. The system was then solvated 
in water resulting in a box of linear dimension 167 A in 
the z-direction. An energy minimization procedure was 
then run, the system was heated up to 295 K, and fi- 
nally the production run started. The first 600 ps were 
discarded to remove artifacts from the initial conditions 
and the information from the subsequent 2 ns collected. 
Other simulation details are as in Ref. [^. 

[75] We calculated the density of water surrounding each ion 
by placing 1 A shells concentric with the ion. The inner 
shells have a larger width to give the same volume. We 
then counted the number of atoms (either hydrogen or 
oxygen) within each of the shells throughout the 2 ns 
simulation at time intervals of 200 fs. Due to the smaller 
bin sizes, the plots have minor differences from Ref. i4L 
at small distances from the ion. 

[76] The electric field was calculated by summing the con- 
tributions from the ion and all partial charges (on the 
hydrogen and oxygen of the water) within 15 A from ev- 
ery field point. The angular component to the field was 
several orders of magnitude smaller because the time- 
averaged field has essentially spherical symmetry. In Ref. 
(41| all water molecules were modeled as dipoles. 

[77] For each time step, all water molecules within a cylindri- 
cal annulus coaxial with the -I- z-axis were taken, where 
the annulus has a 1.5 A width, 1.5 A height, and a central 
ring 5 A from the ion. Then the unit vector connecting 
the oxygen atom of those molecules to the midpoint be- 
tween the hydrogen atoms (the unit dipole pit)) was gen- 
erated together with the unit position vector of the water 
molecules at the centroid of the molecule f{t). We then 
took the scalar product p ■ f per molecule, and averaged 
over the molecules. This set of data was then made into 
a histogram of 501 bins evenly spaced from -1 to 1. 

[78] Here we deal with constant volume and temperature and 



thus use the Helmholtz free energy. 

[79] The most detailed information regarding ion channels 
and physical processes in nanopores is provided by Molec- 
ular Dynamics (MD) - but MD simulations are not able 
to reach the necessary time scales required to extract 
the full information on the current. Indeed, there is a 
hierarchy of approaches going down from macroscopic 
to microscopic models: continuum models - Poisson- 
Boltzmann, Poisson-Nernst-Planck; Brownian dynamics; 
classical then quantum Molecular Dynamics. In practice, 
some combination of the different approaches is often 
used, such as calculating structural/energetic properties 
from MD and using them to construct simpler model sys- 
tems that can then be tested experimentally. This is the 
approach we have followed in this work. 

[80] One may worry that these charge layers - which mainly 
are due to excess ionic density - invalidate the assumption 
of constant ionic density outside the pore. A quick esti- 
mate of the excess density comes from the surface charge 
(for two parallel plates) necessary to give a typical po- 
tential drop of 100 mV over 100 A. This is cr = eoV/l « 
5.5 X 10^'' e/ A . This surface charge is likely contained 
in a layer ~ 10 A thick, giving a density 5.5 x 10~^ e/ A . 
For comparison, at 1 M concentration the density is 
6 X lO"** e/ A - that is, orders of magnitude larger than 
the variation in density necessary to give the electric field 
over the pore. However, increasing the bias or decreasing 
the bulk concentration or pore length may invalidate this 
assumption. 

[81] We note that for all layers to be present, the applied field 
can not be stronger than the ion's field of ~ 0.3 V/A - the 
magnitude of a monovalent ion's field within the third 
layer (~ 7A from the ion) - and approximately double 
that for divalent ions. In this way, the hydration layer 
structure will not be significantly perturbed. 



